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ABSTRACT 



We consider the pair correlation functions of both the order parameter field and its 
square for phase ordering in the 0(n) model with nonconserved order parameter, 
in spatial dimension 2 < d < 3 and spin dimension 1 < n < d. We calculate, in 
the scaling limit, the exact short- distance singularities of these correlation functions 
and compare these predictions to numerical simulations. Our results suggest that 
the scaling hypothesis does not hold for the d = 2 0(2) model. 



1 Introduction 



When systems described by an order parameter with a continuous symmetry are 
quenched from a disordered phase into a broken symmetry phase, the phase order- 
ing dynamics exhibited differs from that of systems with a discrete order parameter 
in a number of important ways. Many different systems in nature possess such sym- 
metries, such as liquid crystals, superfluids etc., and consequently there has been 
considerable interest recently in adapting methods originally developed to investi- 
gate discrete systems to encompass these new symmetries. 

A major goal in the study of systems with both discrete and continuous symme- 
tries, which has been partially fulfilled in recent years, is the approximate determi- 
nation of the scaling function f(x) for the two-point correlation function. Although 
exact calculations have proved possible only for physically uninteresting values of 
d and n |IJ, the spatial and spin dimensionalities, approximate theories [2-7] have 
been proposed which appear to agree closely with simulation results, both for the 
two-point correlation function of the order parameter field 0, 

C 2 (r,t) = (0(x,t)0(x + r,t)> = f(r/L(t)) , (1) 

and its Fourier transform, the structure factor 

S 2 (k,t) = L(t) d g(kL(t)), (2) 

where L(t) is the characteristic length scale at time t after the quench, the angular 
brackets represent an average taken over initial conditions, and the expected scaling 
forms have been anticipated in (1) and (2). In addition, the two-point function of 
the square of the field, 

Cf(r,t) = (0 2 (x,t)0 2 (x + r,i)> - (0 2 (x,t))(0 2 (x + r,t)> , (3) 

(where U denotes 'unnormalized') has recently been calculated within the framework 
of these approximate theories M. For scalar systems, the domain walls are sharp 
so this correlation function is equivalent to the defect-defect correlation function 
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(p(x, £)p(x + r, t)) — (p) 2 , where p(x, t) is the defect density at x at time t. Another 
case of physical interest is that of a system with a complex scalar order parameter, 
e.g. superffuid helium, which can be regarded as an n = 2 vector system. Here the 
correlation function of the order parameter 0(r, t) is not experimentally accessible. 
Physical probes couple to (j) 2 {r,t) instead and so C%(r,t) is the relevant correlation 
function to consider. 

It is generally accepted that, for systems with a non-conserved order parameter 
with n < d, topological defects play an important role in the process of phase 
ordering. For the 0(n) model, one finds that the stable defects are interfaces for n=l, 
vortices for n=2 and monopoles for n=3 ||. The idea we will use in the calculations 
presented in section ^| is that the order parameter field close to one of these defects is 
determined solely by the geometric properties of the defect. Bray and Humayun have 
previously used this observation to calculate the singular short distance behaviour of 
the two-point correlation function, which is responsible for the power law tail in the 
structure factor, and are therefore able to give an exact prediction of the amplitude 
of this tail as well as confirming the result S(k,t) ~ L- n Ar (d+n) @, §. The 



latter generalizes the familiar 'Porod's Law' [TT| to vector systems. We use similar 
ideas to calculate the short distance form of C^r, t), the correlation function of the 
square of the order parameter field, and in section |3.4| we compare simulation results 
to these predictions. We also discuss the validity of the dynamic scaling hypothesis 
in relation for the d = 2 0(2) model, both in relation to our short distance results 
and the simulation data. 

2 Exact results 

In this section we present exact calculations of the short distance behaviour of the 
correlation functions discussed above. In principle, the idea behind the calculations 
can be applied to most correlation functions, although there seems little purpose in 
doing so as in general such correlation functions are not experimentally accessible. 
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The calculation of the singular behaviour of the C^r, t) function at short distance 
has, as we noted in the previous section, already been carried out by Bray and 
Humayun, so we will only summarize the calculation here and refer the reader to 
reference |TU|] for the details. 

Consider the field <fi at points x and x + r in the presence of a point defect at 
the origin, appropriate to the situation n = d. Let us denote the defect density and 
average defect separation at time t by p(t) and L(t) respectively, and the defect core 
size by a. Then for a <C |x| -C L and a <C |x + r| <C L the field (f) at the two points 
will be saturated in length and can be taken to be directed radially outward from 



the origin [0. We therefore have 

7/ \ 7/ \ x • (x + r) 



(4) 
|x||x + r| 

The singular short-distance behavior arises from the rapid spatial variation in the 
direction of the vector near the defect. Holding r fixed and averaging over all 
possible positions of the defect relative to the point x gives 

C S i n g(r,t) = p(t) [ d n x ( j r| X — r ^ — analytic terms] . (5) 

The integral in equation (^) is understood to be over all x which satisfy |x| < L. 
However, the subtraction of 'analytic terms' from the integrand (i.e. analytic in r), 
sufficient to converge the integral at large |x|, enables us to evaluate the integral 
simply by extending it over all space. 

After some algebra the expression for C S i ng obtained is 

C S in g (r, t) = nr n P {t)ir n ' 2 - l B (^y^, du M ~" /2 " 1 { e ' U ~ anal y tic terms ) 

= nry(t)^B (^±1, ^±1) T(-n/2). (6) 

Again, the 'analytic terms' above refer to sufficient terms in the expansion of exp(— u) 
in powers of u to converge the integral. In equation (6), T(x) is the gamma function, 
and B(x,y) the beta function. 

For even n, the gamma function in equation |^ contains poles. To deal with these 
cases, we set n = 2m + e and expand the resulting expression, looking for the order 



unity term in e. We will discuss this procedure more fully below in the context of the 
Ci(r,t) function. Here we simply note that the principle effect is the introduction 
of an additional factor of lnr. 

— * 

To generalize this result to n < d is straightforward. In this situation, (/>(r, t) lies 
in a plane normal to the defect and so we resolve r into two parts; the component 
in this plane, rj_, and the component parallel to the defect core, ry. The average 
over the orientations of is simply the calculation for the point defects above, and 
so the required result comes from averaging equation @ over the orientations of ry, 
with r replaced by its projection onto the plane, rj_. The final result for the singular 
short distance behaviour of the two-point correlation function for general n and d 



is 
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C- fr t) ^ r(-n/2)r(rf/2)I*((n + l)/2) 

Csmg(r ' t) - 7r r((d + n)/2)r(n/2) P{) ' {) 

Once again, for even n we expand the gamma function to extract the leading singular 
term, which leads to a factor of lnr. Taking the Fourier transform of equation (|7|) 
gives the large k form of the structure factor, valid for general n < d, 

S(k , t) = <?M (4^«>/* r2(( " + r ^ 2 2 ' r(<i/2) t-'^» , (8) 

valid for ka <C 1 <C kL(t). We will compare this result to simulation data later. 
Before we do so, let us consider the short distance behaviour of the 2 correlation 
functions. Since 2 saturates to unity at late times, however, it is in practice more 
convenient to consider the pair correlation function of 1 — 2 . In particular, we are 
interested in this function normalized by its large distance limit, 

C 4 (r, f)= <' 1 -^ t »( 1 -?( X + *■'»>. (9) 
(l-^(x,i)>(l-^(x + r,i)) 

In soft spin models of Ising systems (n = 1), the domain wall profile saturates 
exponentially fast to ±1 and so the interfaces are sharp. This makes the calcula- 
tion of C±(r,t) particularly simple. At any point in the system, ((1 — 2 (x, £))) is 
zero unless it is close to a domain wall, so we have immediately ((1 — 2 (x, £))) = 
ap(t) where this relation defines a as the domain wall width. For the product 



(1 — 2 (x, t))(l — 2 (x + r, t)) to be non-zero requires there to be an interface at x 
and at x + r. If r is small compared to the characteristic distance between interfaces, 
L(t), then the two points must lie on the same interface. If we consider the point at 
x to lie on an interface, then the additional requirement that x + r lies on the same 
interface gives a factor aSd-i/ Sjr, where Sd is the surface area of a d dimensional 
unit sphere. We have, therefore, for n = 1, 

((1 - 2 (x, f))(l - 2 (x + r, t))) = a 2 ^^ , a « r « L(t) , (10) 

which gives, for n = 1 and general d, 



The d = 3 result was recently derived by Onuki [|13 

For vector systems, the field far from a defect saturates as 1 — 2 (x) = a 2 /x 2 , 
where x is the distance from the defect, and a is the 'size' of the defect core. It 
follows, therefore, that the short distance behaviour of C 4 (r,t) is given by 

<(l-0 2 (x,t))(l-0 2 (x + r,t))> 



C 4 (r,t) 



((l-0 2 (x,t)))((l-^ (x + r;t))) 

f^/ Lrfx x^)/^/ Lrf 4) 2 ' ^ 



provided the integrals are dominated by the field close to a defect. As before, we 
consider first the case of point defects, n = d, then generalize the result to n < d. 



The denominator in equation (|12|) can be evaluated by changing to polar coordinates 
and gives, for n = 2, 

Cf = (2vra 2 pln(L/a)) 2 , n = 2 . (13) 

For n > 2 the denominator in (|1^) is no longer dominated by the field close to 
the defects and so we cannot treat these cases in the same detail as for n = 2. 
However, we can still evaluate the leading short-distance behavior of the numerator 
for 2 < n < 4, and the leading short-distance singularity of the numerator for any 
n > 2. 



To evaluate the numerator (denoted by superscript W), we first make a rescaling 
x — > rx which gives 

Cf(r,t) = ((l-<?(x,t))(l-0 2 (x + r,t))> 



a 4 



= pr™- 4 / cix (14) 

Ja/r<x<L/r \X.\ Z \X. + T\ z 

We now take L/r — ► oo and a/r — > 0. For n = 2, this makes the integrals become 
logarithmically divergent, and the following treatment applies strictly only for 2 < 
n < 4, as we shall see below. We will discuss how the calculation is extended to 
n = 2 later. 

Using the integral representation J °° dv exp(-uv) = 1/u gives 

roc roc roo 

Cf (r,t) = pa 4 r"" 4 / du / dv / dx exp(-w|x| 2 - v\x + f| 2 ) . (15) 

JO JO J-oo 

Evaluating the x-integral first, followed by the (u,v) integrals, gives, after some 
algebra, 

= ^ r((4-yi(,-2 W 

1 (n — 2) 

To extend this result to n < d we consider the vector r to be resolved into two 
components, and ry, where is the (d — n)-dimensional component of r in the 
plane perpendicular to the defect core and rn is the orthogonal vector containing 
the remaining components of r. Using generalized spherical polar co-ordinates, the 
projection of r into the plane perpendicular to the defect is simply 

d—n 

|r_i_| = r Yl I shift; | (17) 
i=i 

To calculate C^(r, t) we use the result from equation (|i"6|) with r replaced by and 
average over all orientations of r, which is equivalent to averaging over the possible 
relative positions of the two points to the defect. This gives 

1 [n — 2) 

To evaluate this expression we require the average 

Id- n \ 
|rx| n ~ 4 ) = r n ~ 4 (J] I sin ftp- 4 ) 



ii=l 
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d-n / n / 2 / r 7r/2 N 

" H ' A II / d0i(sin0 i ) ft ~ 4+<i ~ 1 ~ i / / ^(sin^)^ 1 -' 



i=i 



o / Jo 



= 4 r(n-2)r(d/2) 

r((n + d)/2-2)r(n/2)" 1 ; 

Inserting this into gives, for general d > n and 2 < n < 4, the short- distance 
behavior 

C N (r t) - T n/ 2 r(2-n/2)r(d/2)r^(n/2-l) ^ 

C4(r ' t)_7r r((n + d)/2-2)r(n/2) ' a « r « L W • (20) 

In fact this result gives the leading short-distance singularity in for all n > 2 
(with the usual logarithmic factors for even n, which can be extracted by taking the 
limit from real n appropriately). For 2 < n < 4, it gives the leading short distance 
behavior (not just the singularity), whereas, for n > 4, O4 saturates to a constant 
at short distance. (By 'short distance' we mean, of course, a <r <C L(t) as usual). 
Equation (20) implies the power-law tail S^k, t) ~ p(t) / c -( a! + n - 4 ) ) for d > n > 2, in 
the Fourier transform of C4, a form first suggested in ||. 

For n = 2, there is a double pole in equation fl2"0| ) for d 7^ 2 and a single 
pole if = 2. These poles come from the lower end of the integration range in 
equation flTiD . In reality, (1 — 2 (x, t)) = a 2 /x 2 breaks down at x ~ a so the 
integrals in equation ( |T4"D should be cut off at r ~ a. The effect of introducing the 
cutoffs would be to eliminate the poles in equation (0). To extend this result to 
n = 2 then, we set n = 2 + e and expand the resulting expression to find the 0(1) 
term in e. The terms which are divergent as e — ► would be removed by a short 
distance cutoff. This gives, for d > 2, 

#f*« = r(i- e /2)rK/2)ry 2 ) 



Cf(r,t) = p(t)a 4 vr 



r((e + d)/2-l)r(l + e/ 2 )r 2 - 
(d - 2)7rp(t)a 4 ln 2 (r/a) 



(21) 



on retaining the 0(1) term in e. We introduce factors of the form ln(r/a), rather 
than ln(r/L) say, because the cutoff required on the integrals is at r = a. For d = 2, 
we need to take the limit d — > 2 before taking n — > 2. This leads to a cancellation 
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of two of the gamma functions in equation (pOf ) leaving 

27rp(t)a 4 ln(r/a) ^ 



on retaining the 0(1) term in e. A direct calculation of the n = 2 case with cutoffs 
on the integrals confirms these results, which are correct to leading logarithmic 
accuracy. 



In the scaling limit, where r/L is kept fixed and L — > oo, equations (|21f) and 
( p2|) can be simplified using ln(r/a) = ln(r/L) + ln(L/a) ~ ln(L/a). The normalized 
correlator O4 = C^/Cf is therefore, using fll3|) for Of, 



tZ-2 , , . 

= 4^' rf>2 ' < 24 > 

valid for a r L(t). It is interesting that for d > 2, the logarithms in and 

Cf cancel to produce an expression in a conventional scaling form (with scaling 

variable pr 2 ), as would be expected, whereas for d = 2 a single logarithm remains 

in the denominator and therefore C±(r, t) does not have a scaling form at short 

distance ||. This does not necessarily imply that there are no scaling phenomena 

in the d — 2 0(2) model but it does show that at least this correlation function 

does not scale. Taken together with evidence from the simulations of |f[4 |, and 



our own simulation results presented below, that indicate a multiplicity of length 
scales in this system, as well as the calculation of Bray and Rutenberg of the energy 



dissipation [15], it now seems likely that scaling does not occur in the d = 2 0(2) 
model. 

In summary, we have exact predictions for the short distance singularity of 
02(1", t), or equivalently, for the tail amplitude of the structure factor. We have 
also predicted the short distance form of O^r, t) for n < 2, which for the 2d XY 
model violates the scaling hypothesis. We now compare these results to computer 
simulations. 
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3 Simulation results 



We now turn our attention to computer simulations of the 0(n) model, the results 
of which we will compare to the above predictions. For n > d + 1 there are no stable 
topological defects and for n = d + 1 the defects are extended textures, the dynamics 
of which are not well understood. We therefore have restricted our investigations 
to 2 < d < 3 and 1 < n < d. The data from the n = 1 simulations, some of which 



have already been published in another form jig] , were provided by S. Sattler. 
3.1 Cell dynamical simulations 

The simulation of phase ordering using cell dynamical simulations (CDS) is by now 
standard ||17|| . The essential idea is to take a lattice of 'soft' spins, corresponding to 
a coarse-grained order parameter field <j)(r,t), and update them with the mapping 

$ n+1 (t) = D ^^Mj) - +E$ n (i)tanh(\$ n {i)\) (25) 

where = (f>/\(f>\ (for scalar fields <fi = sgn(0)), z is the number of nearest neighbours 
and D and E are adjustable parameters, for which we chose the standard values 
D = 0.5 and E = 1.3 |14]]. One can regard equation (|2~5D as a discrete TDGL 
equation in which the first term is the discretized Laplacian and the second is derived 
from a potential V(4>), the form of which has been chosen to keep the iteration 
process numerically stable and to give the same attractors of the dynamics as the 
TDGL equation. 

It is desirable when measuring C2 to 'harden' the spins, which is achieved by 
computing (for general n) 

C?(r,t) = ($(x,t)4(x + r,t)). (26) 

This reduces the effective size of the defects (domain walls, strings etc.) to zero 
and hence scaling can be achieved at smaller values of r. Unfortunately, the same 
procedure cannot be applied to C4 (which would vanish identically) and hence the 
scaling in this case for small r is considerably worse. 
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To measure the defect density in the vector systems we used two different algo- 
rithms. For the d = n systems, the defects are points and so can be identified by 
looking for peaks in the local energy density = — J2<j><fii • <t>j, where the sum 
is over the nearest neighbors of i. For the d = 3 0(2) model, however, the defects 
are strings and consequently this method is not so simple to use. Instead, we look 
for plaquettes around which the order parameter rotates through ±2ir radians and 
associate with each one a unit length of string. The same approach applied to the 
d = 2 0(2) model gives results identical to those obtained with the algorithm based 
on the energy density. It is important to note, however, that the length of string 
given by this procedure is that given by a Manhattan metric, because the string is 
constrained to lie on the dual lattice. To compare to our calculations in section [2|, 
which were based on a continuum description, requires that we use a Euclidean met- 
ric to measure the string length. Assuming that on a scale large compared to the 
lattice spacing the strings are isotropically oriented, we can convert between these 
two lengths by adjusting by a factor [Lm/^e] = 3/2 for d = 3, where the brackets 
represent an average taken over all orientations of the string and Lu and are the 
Manhattan and Euclidean lengths of the string. 

The details of the simulations are shown in table In all the following plots, we 
will use the symbols circles, squares, triangles, diamonds and crosses to represent 
progressively later measurement times. The actual times to which these symbols 
refer can be read off from the table. 

3.2 Dynamic scaling 

We start the analysis of the data by considering a test of the dynamic scaling hy- 
pothesis. If this hypothesis is true, then the data for the two-point correlation 
function C^r, £) can be collapsed onto a single curve for each system by rescaling r 
by the characteristic scale L(t). This scale can be defined in a number of ways. For 
example, it could be defined as the value of r for which C^r, t) = 1/2, or even as 
the value which gives the best overall collapse of the data. A definition of this sort 
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Table 1: Details of lattice size and number and length of runs for the numerical 
simulations. 



can be rather misleading, however, because a free parameter has been introduced for 
each value of t. As the form of C2(r, t) is rather featureless, it is not too surprising 
that a good collapse of the data can be obtained using methods such as these. An 
alternative and superior method is to use some characteristic length scale that is not 
directly related to C 2 (r, t). If the scaling hypothesis holds then the defect density 
should scale as pit) ~ L~ n (t), and so it seems logical to use rp 1//n (t) as the scaling 
variable. For the n = 1 systems we have used, instead of p, (1 — 2 ) which scales in 
the same way as the true defect density. Figures [l] to || show the two-point corre- 
lation function scaled in this manner along with the theoretical predictions of OJK 
for the scalar systems and of Bray, Puri and Toyoki (BPT) for the vector systems. 
The general form of these results is [0, [| 

'71+ 1 l^ 2 



Ztt 



B 



'.2'2' 2 ' 7 



(27) 



2 ' 2, 

where F(a,b; c; z) is the hypergeometric function, and 7 = exp(— x 2 /8), with x = 
r/L(t) the scaling variable. In the theoretical curves, the horizontal scale has been 
adjusted to give the best fit. 
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Apart from d — 2, n — 2, the collapse of the data is excellent. For d — 3, n — 3 
there is some evidence for finite size effects in the last data set, but otherwise the 
data scale well. We conclude that our numerical results are entirely consistent with 
the scaling hypothesis in these systems. It is clear, however, that for the d = 2 
0(2) model, i.e. the 2d XY model, the data do not scale. Furthermore, there is no 
evidence that the data at later times are approaching an asymptotic scaling curve. 
It is possible that the system has yet to reach the scaling regime even at the latest 
times, but considering how early the scaling regime is reached for the other systems, 
it is also possible that the scaling hypothesis fails for this system. It is interesting 
that the data can be collapsed to a single curve by using values for L(t) specially 
chosen to do so. The quality of the scaling achieved is as good as in any of the other 
systems (see figure D). Bray and Rutenberg also find evidence that there may be 
scaling violations in the 2d XY model [|15|]. They conclude that energy dissipation 
due to the ordering process occurs significantly on all scales between the vortex 
core size and the inter-vortex spacing, suggesting there may not be a single length 
scale L(t) that characterizes the system morphology. We recall that the exact short 
distance behaviour of C^r, t) calculated in the previous section is also not in a 
scaling form. In all cases, the theoretical curves fit very well with the data. There 
are small discrepancies, but considering the approximations made in the calculation 
of these curves, it is rather surprising how close they appear to be to the simulation 
data. We will see later, when we eliminate the free parameter required to compare 
the theory to simulation results, that there are in reality considerable discrepancies 
which are hidden by the introduction of the free parameter. 

3.3 Growth laws 

The growth law for the characteristic scale in systems with a non-conserved order 
parameter is expected to be L(t) ~ t 1 ^, with z = 2 []18 , [15[| . By fitting a straight 



line to a log-log plot of the defect density against time we find the values of l/z 
shown in table H which follow from the relation L(t) ~ p(t)~ n . After infinite time 
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0.46 
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Table 2: Growth laws for the characteristic scale in 0(n) systems. The values were 
obtained from the decay of the defect density. 

all the defects should have been removed from the system. Consequently, we expect 
p(t) plotted against t~ n l z to be a straight line passing through the origin. Figure [7] 
shows these plots along with a straight line fit. To quote errors on the values of 
1/z is somewhat misleading because the data at different times come from the same 
system and are therefore highly correlated. In all cases, however, a fit with z = 2 is 
unacceptable, even for the n = 1 systems! The reason for this is unclear. It could 
be that the systems have yet to reach their asymptotic behaviour. A more likely 
explanation is the weak pinning effect on the defects, introduced by the lattice, 
which becomes more important at later times as the driving forces become weaker. 
However, the data are consistent with zero defect density for t — > oo, so there is 
no direct evidence that defects become frozen at late times, as they do for hard- 



spin vector simulations |19| . It should be noted that previous simulations on vector 
systems also find growth exponents consistently less then 1/2 [14], |20[| . 

For the 2d XY model, the growth laws for p(t) -1 / 2 and L(t), the length scale 
required to collapse the C*2(r,i) function, are noticeably different. The log-log plots 
shown in figure |8| show clearly the different power laws involved. A least squares fit 
gives p(t) -1 / 2 ~ t ' 37 and Lit) ~ t 0A2 . An estimation of the errors on these expo- 
nents is difficult due to correlations between data points. Mondello and Goldenfeld 
have previously simulated the 2d XY model with rather longer run times than we 
used ||14j| . Their data also provides evidence for different growth laws for different 
length scales for the time ranges they used. 



13 



3.4 Short distance results 

3.4.1 The structure function tail amplitude 

We now consider the predictions of Bray and Humayun for the amplitude of the 
power law tail in the structure factor. We present only our results for vector systems 
as a comparison to Ising simulations has already been made |I(J. The scaling form 
for the structure factor is S(\s.,t) = L d (t)g(kL(t)) and we expect a power law tail 
of the form L[t)~ n k~^ d+n \ For large k, S(k, t)k d+n / p{t) should therefore be the 
constant predicted by equation (|]). Using L{t) = p~ 1 l n (t), the scaling variable 
is kp~ l l n (i) which we plot against S(k,t)k d+n / pit) = h(kp~ 1 l n (t)), where h(x) is a 
scaling function. Figures |9| to [11] show these plots along with the analytic predictions 



for the large /cp _1 / n limit. Note that we have plotted (S(h,t)k d+n / p(t)) rather than 
(S(k, t)k d+n ) I (p(t)) because the calculated amplitude of the tail should be correct in 
each individual system. In all three cases, the numerical data is consistent with the 
predicted values of the tail amplitudes. For the largest values of k, there is a tendency 
of the data points to move up away from the predicted line. This is presumably 
because the calculations were carried out using a continuous order parameter field 
and at these large k, the effect of the lattice becomes noticeable. The case d = 3, 
n = 2 shows a small but systematic discrepancy from the predicted value at large k. 
Again we attribute this to a lattice effect. In particular, the correction applied to 
the data to compare lengths measured on the lattice (Manhattan metric) to those 
calculated on the continuum (Euclidean metric) becomes inappropriate at lengths 
scales comparable to the lattice spacing. This effect accounts for the sign of the 



discrepancy in Figure 10 



An interesting feature of the simulation curves is that they always approach their 
asymptotic constant value from above. For a scalar order parameter, Tomita |21 



has accounted for this in terms of the curvature of the interfaces, but we are not 
aware of any similar calculations for vector fields. This same feature can also be seen 
when the data is plotted in the more usual form \n(L~ d (t)S(k, t)) against ln(fcL(t)), 



14 



although, because of the logarithmic scale, the effect is less pronounced. For com- 
pleteness, we present the data replotted in this form in figures [12] to [14], along with 
the theoretical curve predicted by BPT. The rescaling factor introduced into the 
theoretical curves in these plots was the same as that used to give the best fit to 
the two-point correlation function in real space. Because this factor was determined 
using the whole of the two-point correlation function, the curves fit the data best 
at small k. Note that, since the data approach their asymptotic straight lines from 
above, there is a tendency in fitting data to assign too negative a slope to the 
structure factor tail. 

3.4.2 C 4 for the scalar systems 

For n — 1 and n = 2 we have exact predictions for the short distance behaviour 
of C4(r,t), the correlation function of the square of the order parameter. Plotting 
C<i(r, t) against C^r, t) gives a parameter-free scaling curve that can be used to test 
theories of phase ordering JTB[. Here we will consider the short distance behaviour 
of these functions. From the previous section, we have for n = 1 

C 2 (r,t) = 1 - 2T ~' / V((d+ / 1V2) P{t)r ' a « r « L W- 

^ = wr ■ < 28 > 

and so 

^(M)^ r((rf -™ +1)/2) (l-^(r, t )) (29) 

in this regime. Therefore we expect a plot of C 4 T 1 (r,t) against 1 — C2(r,t) to be 
linear near the origin with gradients 7r 2 /4 and 2 for d = 2 and d = 3 respectively. 



Figures [15] and 116] show the data plotted in this manner for the 2d and 3d n = 1 



systems, along with a line of the appropriate gradient [22|. The data depart from 
the expected line for very small r because we require a<r< Lit). It does appear 
that for asymptotically large times the data will lie on the line we predict. It is also 
possible to calculate C^ir.t) as a function of 1 — C 2 (r,t) using the approximate 
theories of OJK 0] and Mazenko || mentioned above. From equation (27) with 
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7? 



1 one obtains C 2 (r,t) = (2/tt) sin" 1 ^) while g |16| C 4 (r,t) = 



where 7 is the normalized two-point function for a gaussian auxiliary field (see 



references || and [I6|). Eliminating 7 gives the parameter-free relation 



7T 



C^ = sin -(1-C 2 ) . (30) 



This prediction is included in Figures [15] and [16]. Plotting the data this way provides 
an absolute test of theory. Clearly, the existing theories based on a gaussian auxiliary 
field are not quantitatively accurate. It is interesting to note, however, that for small 



r, i.e. small I — C2, equation (|30|) gives the same linear relationship between C 4 and 



I—C2 as the exact small r result in equation ( p9|) does, and, moreover, the coefficients 



are the same in the limit d — > 00. It is also clear from figures [15] and that the data 
are closer to the predicted form in equation (|30"D for d = 3 than for d = 2. These 
facts add to the growing body of evidence that the OJK and Mazenko theories are 
exact in the large d limit p3| . Bray and Humayun |23| have recently discussed a 
procedure by which the leading order OJK result can in principle be systematically 
improved. The present data will provide a good test of any improvements to the 
theory. 

3.4.3 C 4 for the XY model 

From section ^] we have, for n=2, the results given by equations ( fZ3"| ) and ( |2~3j ) for 
C^(r,t) for d = 2 and d > 2 respectively. Figure [17| shows the simulation data for 
C^r.t) for the d = 3 XY model plotted against the expected scaling variable pr 2 . 
The data scale well for large r and are consistent with the predicted short distance 
behaviour, C 4 ~ 1 (r,7;) = Anpr 2 . 

We noted above that we cannot calculate C^r,/;) for the 0(3) system because 
(1 — (f) 2 (r,t)) is not determined by the field close to the defects. The result (|20"D for 
C^(r,t) = ((1 — 2 (x,i))(l — 2 (x + r,r;))) is still valid, however, although we cannot 
use it as an absolute prediction because the defect core size a is not known. From 
equation ([H]) we expect, for n = 3, C4 ~ 1/r for small r. The scaling function for 
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1/C 4 should therefore be linear for small r, which is certainly consistent with the 
data shown in figure [18] (recalling that the data break away from the true scaling 
curve for small r, when the condition r>ano longer holds). 

For the d = 2 XY model, we predicted in section |2| that C±(r, t) would not scale 



for small r, but instead would have the form given by equation (|23|). Figures [L9 



and |20] show the data plotted in the form predicted, and against the naive scaling 
variable pr 2 . In both cases the data appear to scale reasonably well. There is an 
ambiguity in the definition of L(t), the characteristic scale, used in figure [l^, so we 
have chosen L(t) = l/y/2p, simply because this gives the best agreement with the 
prediction for the small r behaviour. Different choices of L(t) still give C^ 1 oc r 2 for 
small r, but with different gradients. Our simulation results are therefore consistent 
with the predicted form of C±(r,t), although they do not constitute an absolute 
test. It is interesting that, in contrast to the two-point function, the data can 
be scaled using V ' pr' 1 . The deviations from scaling predicted by equation ( p3|) are 
only logarithmic in L(t), however, which does not vary greatly for the range of L 
available from the simulation. The fact that C^(r,t) can be approximately scaled 
using the defect density, whereas C 2 (r,t) cannot, may be because C 4 (r,t) is more 
closely associated with the defects. For n — 1 systems, for instance, because the 
domain walls are sharp, C 4 is essentially a defect-defect correlation function. 

Just as for the scalar systems, we can plot 1 — C^r, t) against l/C^r, t) to obtain 
a parameter-free test of the BPT result |2|, [|, equation (|27D , and the related theory 
of references [0] and 0. The equivalent result for C 4 was obtained in reference [8]: 
C 4 = F(l, 1; n/2; 7 2 ). Eliminating 7 between C 4 and C 2 gives C^ 1 as a function of 
1 — C*2. Figures [TL] to |23| show the data plotted in this form [22]. The data do not 
fit at all well quantitatively to the theoretical prediction, but do improve for larger 
d, as expected [|3[ . 
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4 Conclusions 



We have performed numerical simulations of the 0(n) model for 1 < n < d in 
dimensions d = 2 and d = 3. Exact calculations for the short distance behaviour of 
C*2 and C4 have been compared to the simulation data and satisfactory agreement 
found. A particularly interesting feature of these exact results is that for the 2d XY 
model they predict a form for C4 which cannot be expressed in a scaling form, a fact 
already noted by Bray in the context of his approximate calculation of C4 |§. Also, 
C2 does not scale when plotted against rp(t) 1//2 . We therefore suggest that there are 
logarithmic violations of scaling in this model. 
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Figure captions 

1. The two-point correlation function of the d — 2 scalar model. The data have 
been scaled using (1 — <p 2 ) which is proportional to the defect density (data 
courtesy of S. Sattler). The solid curve is the OJK prediction, with the hori- 
zontal axis scaled to give the best fit by eye to the data. 

2. The two-point correlation function of the d = 2 0(2) model. The data have 
been scaled using the defect density. 

3. Same as Figure 1, but for the d = 3 scalar model. 

4. The two-point correlation function of the d — 3 0(2) model. The data have 
been scaled using the defect density. The solid curve is the BPT prediction, 
scaled to give the best fit by eye. 

5. Same as Figure 4, but for the d — 3 0(3) model. 

6. The two-point correlation function for the d — 2 0(2) model. The data have 
been scaled using the values L(t) which gave the best scaling. The solid curve 
is the BPT prediction, scaled to give the best fit by eye. 

7. Growth laws for the defect density in 0(n) models. Moving left to right and 
top to bottom, (d,n) takes the values (2,1), (3,1), (2,2), (3,2) and (3,3). 

8. Log-log plot of p(i) _1//2 (top) and L(t) (bottom) against time for the d = 2 
0(2) model. The straight lines have gradients 0.37 and 0.42 (top to bottom). 

9. Structure function of the d = 2 0(2) model, with the power-law tail scaled 
out. The line is the exact result of Bray and Humayun for the amplitude of 
tail. 

10. Same as Figure 9, but with d = 3, n = 2. 

11. Same as Figure 9, but with d — 3, n — 3. 
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12. Structure function of the d = 2 0(2) model. The line is the curve calculated 
by BPT, scaled to give the best fit by eye. 

13. Same as Figure 12, but with d = 3, n = 2. 

14. Same as Figure 12, but with d — 3, n = 3. 

15. Plot relating the correlation functions ^(r, t) and C±(v,t) for the 2d scalar 
model (data courtesy of S. Sattler). The straight line is the exact small r 



prediction and the curve is the 'OJK' prediction || JL6 . 

16. Same as Figure 15, but for the d — 3 scalar model. 

17. The correlation function d(r, t) for the 3d XY model. The straight line is the 
predicted small r behaviour. 

18. The correlation function C^(r,t) for the d = 3 0(3) model. The small r 
behaviour is expected to be linear at asymptotically late times. 

19. The correlation function C^(r, t) for the 2d XY model. The straight line is the 
predicted small r behaviour. We have used L(t) = k^/p where the parameter 
k was chosen to give the required gradient. 

20. The correlation function C±(r, t) for the 2d XY model plotted in the naive 
scaling form l/C^r, t) against pr 2 . 

21. Plot relating the correlation functions C^r, t) and C±{Y,t) for the 2d XY 
model. The curve is the 'BPT' prediction || [16|]. 

22. Same as Figure 21, but with d = 3, n = 2. 

23. Same as Figure 21, but with d — 3, n — 3. 
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